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ABSTRACT 



Context. The origin of solar transition region redshifts is not completely understood. Current research is addressing this issue by 
investigating three-dimensional magneto-hydrodynamic models that extend from the photosphere to the corona. 
Aims. By studying the average properties of emission line profiles synthesized from the simulation runs and comparing them to 
observations with present-day instrumentation, we investigate the origin of mass flows in the solar transition region and corona. 
Methods. Doppler shifts were determined from the emission line profiles of various extreme-ultraviolet emission lines formed in the 
range of T = 10* - 10*" K. Plasma velocities and mass flows were investigated for their contribution to the observed Doppler shifts in 
the model. In particular, the temporal evolution of plasma flows along the magnetic field lines was analyzed. 

Results. Comparing observed vs. modeled Doppler shifts shows a good correlation in the temperature range log(r/[K])=4.5-5.7, 
which is the basis of our search for the origin of the line shifts. The vertical velocity obtained when weighting the velocity by the 
density squared is shown to be almost identical to the corresponding Doppler shift. Therefore, a direct comparison between Doppler 
shifts and the model parameters is allowed. A simple interpretation of Doppler shifts in terms of mass flux leads to overestimating 
the mass flux. Upflows in the model appear in the form of cool pockets of gas that heat up slowly as they rise. Their low temperature 
means that these pockets are not observed as blueshifts in the transition region and coronal lines. For a set of magnetic field lines, 
two different flow phases could be identified. The coronal part of the field line is intermittently connected to subjacent layers of either 
strong or weak heating, leading either to mass flows into the loop (observed as a blueshift) or to the draining of the loop (observed as 
a redshift). 

Key words. Stars: coronae — Sun: corona — Sun: transition region — Sun: UV radiation — Techniques: spectroscopic 



1. Introduction 

Persistent average redshifts of em ission lines from the so lar tran- 
sition region were observed by iDoschek et al.l (1 19761) for the 
first t ime. Early observations from rockets (e.g., Hassler et al., 
1199 lb had limited potential and did not provide a conclusive pic- 
ture of the nature of these redshifts. The evaluation of SUMER 
(Sola r Ultra violet Measurements of Emitted Radiation, Wilhelm 
et al.. ll995h observations onboard SOHO led to the conclusion 
that redshifts in the transition region predominately present ver- 
tical motion s that turn into bluesh ifts at higher temperatures 
(Peter, 119991: Peter & Judge, 1 19991) . These redshifts have been 
explained by a variety of models including the re turn o f previ- 
ously heated spicular material (Athay & Holzer, Il982l Athay, 
|1984|) or nanoflare-induced waves generated in the corona that 
would produce a net redshift in transition region lines (Hansteen, 
Il993h . For a more comprehensive discussion of coronal loops 
in rela tion to the D oppler shift o bservations see iPeter & Judge! 
(119991) . Peter et al. (|2004 l2006l) found that three-dimensional 
numerical models spanning the region from the photosphere to 
the corona are producing redshifts in transition region lines. 
In these models, coronal heating comes from the Joule dissi- 
pation of currents produced by the braiding of the magnetic 
field through photospheric plasma motions; however, details 
on the origin of the redshifts have not been identified so far. 
ISpadaro et alJ (l2006h used one-dimensional coronal loop mod- 
els to show that the redshifts can be t he result of transient h eat- 
ing near the loop footpoints. Recentlv. lHansteen et alJ (1201 Oh has 



shown that transition region redshifts are naturally produced in 
episodically heated models where the average volumetric heat- 
ing scale height lies between 50 and 200 km. The authors suggest 
that chromospheric material is heated in place to coronal tem- 
peratures or that rapidly heated coronal material leads to a con- 
duction front heating the chromospheric material below. Both 
options lead to a high-pressure plug of material at upper tran- 
sition region temperatures that relaxes towards equilibrium by 
expelling material. 

In this paper, we describe the investigation of a small ac- 
tive region by means of a three-dimensional (3D) magneto- 
hydrodynamic (MHD) model. We synthesize transition region 
lines that have been typically observed by SUMER. The synthe- 
sized emission may be compared directly to observations. The 
mass balance in the model simulation is investigated as is the re- 
lation of the Doppler shifts to the magnetic structure in order to 
elucidate the nature of the transition region redshifts. 

The 3D MHD model is introduced in section |2] Section [3] 
focuses on the spectral line synthesis from the simulation results 
and on the comparison with observations. The origin and nature 
of the modeled Doppler shifts is investigated in section|4] before 
results are discussed in section |5] 

2. The coronal 3D MHD model 

We are applying a 3D MHD model of the solar atmosphere 
above a small active region in a computational box of the 
size 50x50x30 Mm^. A detailed description of the model can 
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be found in iBingert & Peter! (1201 The simulations are car - 
ried out using the Pencil-code dBrandenburg & Doblerl l2002l) . 
a high-order finite-difference code for compressible MHD. The 
set of equations being solved includes the continuity equa- 
tion, the induction equation, the Navier-Stokes equation, and 
the e nergy equation with an anisotropic Spitzer heat conduction 
term ('SDitzer"l962') and an optical thin radiative loss function 
JCook et al. 1989). 

For the lower boundary we constructed a map of the mag- 
netic field based on two observed magnetograms, one from an 
active region, spatially scaled down by a factor of five, and an- 
other one from a quiet Sun area. Thus, the resulting magnetic 
map represents the quiet Sun with some activity added by a pore 
and is supposed to be typical for the magnetic field underlying 
the transition region of the Sun outside active regions. 

The main difference betwee n the model pre s ented in this 
study and the one described by iBingert & Peter! (1201 ih is the 
higher grid point resolution of our simulation, allowing for a 
slightly lower magnetic diffusivity. The resolution of 256^ grid 
points corresponds to a grid spacing of 195 km in the hori- 
zontal direction and 75-285 km in the vertical direction, where 
a nonequidistant grid is applied. This corresponds to a maxi- 
mum resolvable vertical temperature gradient of about 2 K/m, 
when assuming a temperature difference of 10^ K across five grid 
po ints. This number is confirmed by the 3D model discussed 
by iPeter et aP (l2006l) . see their Fig. 5. In a stat ic ID ernpirica l 
model based on an emission measure inversion, iMariskal (119921) 
found a maximum temperature gradient of the order of 10 K/m 
(his Table 6.1). Flows present in our 3D model are assumed to 
smooth out steep temperature gradients, resulting in temperature 
changes that are less substantial than ex pected for the cas e of the 
static ID empirical models (see Fig. 5 of lPeter et al.l2006l) . Thus, 
it can be assumed that temperature gradients in the 3D model 
are being properly resolved. However, to synthesize transition 
region emission lines, the temperatures and densities as com- 
puted in the 3 D MHD model h ave to be interpolated (Sect. l3.11 
for details see lPeter et"ani2006h . 

The visco sity and heat conduction terms as published by 
ISpitzeil (Il962h are applied; i.e., in the coronal part of the box, 
the dynamic viscosity is set to vx!lO"m^s"' and the thermal 
diffusivity is about ;i'ailO'^ m^s"'. For velocities of up to 100 
km/s and a grid scale as used in our simulation, the average 
Reynolds number Re and the average Prandtl number Pr are of 
the order of 0. 1 . The magnetic diffusivity is chosen in a way that 
the average magnetic Reynolds number Rm is near unity, i.e., 
TyxilO'^m^s"'. This ensures that the current sheets that are being 
formed (i.e., the gradients in magnetic field) can still be resolved 
in the computation. This choice of parameters implies that the 
average magnetic Prandtl number is Pm-v/Ti^lQ, i.e., on aver- 
age, magnetic field gradients are larger than velocity field gradi- 
ents. The above-mentioned average dimensionless numbers can 
only be used for a rough characterization, because they vary con- 
siderably by orders of magnitude in the computational domain, 
depending on the local plasma conditions. 

While the average magnetic Prandtl number is greater than 
unity, implying that the velocity field is more effective in am- 
plifying field on the smallest scales than the magnetic field, still 
the dissipation of energy in the corona is dominated by Ohmic 
heating processes, and viscous heating only plays a minor role. 
The deposited energy is efficiently redistributed by heat conduc- 
tion, the rate of which is comparable to the Ohmic heating rate, 
except for low temperatures, where radiative losses dominate. 

The boundary conditions are chosen such that the simulation 
box is periodic in (x,y) and closed for mass flux on both the top 
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Fig. 1. Height dependence of the average temperature in hori- 
zontal layers of the simulation box (solid line). Diamonds in- 
dicate the formation temperature of the emission lines that are 
investigated in this study (compare Table [T] column labeled 
log Tconirib)- The dashcd lines show the minimum and maximum 
temperature in the horizontal layers. 



and bottom. At the bottom boundary, the temperatures and den- 
sities are set to the average photospheric temperature and density 
values at T=l. The magnetic field at the lower boundary is driven 
by horizontal motions mimicking the photospheric horizontal 
granular motions. This leads to the braiding of the magnetic field 
lines and the induction of cuiTents, which are subsequently dis- 
sipated by Ohmic heating. As a result, the plasma is heated and 
driven, leading to a dynamic corona. The simulation was run - 
for some 30 minutes - to a state where the behavior is indepen- 
dent of the initial condition. This point of time is defined as f = 
min in our study. The simulation produces a dynamic transition 
region and corona with flows up to and above the sound speed 
and with the horizontally averaged quantities, such as the aver- 
age temperature (see Fig.[T} and density as a function of height, 
showing only a little variation over time. Below about 4 Mm, 
the computational domain contains material at and below chro- 
mospheric temperatures. This height encompasses ten chromo- 
spheric scale heights, i.e., a drop in density by more than four 
orders of magnitude. This dense part of the box acts as a mass 
reservoir for evaporating plasma into the corona and for material 
draining from above. It also serves as a reservoir for energy that 
is conducted back from the corona without being radiated in the 
transition region. 

The energy deposition rate in the corona is highly variable 
in space and time. A reasonable estimate of (2~10"'' Wm""* for 
the energy deposition at a height of some 10 Mm corresponds 
to an energ y flux density of some 100 Wm"^ at the base of the 
corona (see Bingert & Peterll201 ll their Fig. 5). A typical length 
scale on which current sheets are forming is 1 Mm. Thus, the 
associated time scale for magnetic dissipation is of the order 
of T=^^/77^100s. Therefore, the deposited energy density can 
be estimated to be E-Qtk\Q'^ ] voT^ . According to the mag- 
netic energy density B^/{2fi()), this corresponds to a field strength 
of about Bx5 G. This number is within an order of magnitude 
of the actual magnetic field strength in the coronal part of the 
computational domain, which is about lOG at 10 Mm height. 
Consequently, the field line braiding in our model can lead to an 
accumulation of magnetic energy within an order of magnitude 
in energy to the corresponding potential field extrapolated from 
the photosphere. If the driving in the photosphere is switched 
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off, the heating will only continue for the duration of the dissi- 
pation time T, i.e., some minutes. Thus, there is a need for a con- 
tinuous energy supply through the photospheric motions. The 
above-mentioned order of magnitude estimations are confirmed 
by a more detailed analysis of the numerical experiment. 

3. Results 

3. 1 . Spectral line synthesis 

Optically thin transition region and coronal emission lines are 
synthesized from the results of the MHD calculations (Table 
[1). In the low-density plasma of the solar transition region and 
corona, spontaneous radiative decay is the main process that de- 
populates excited levels. The rate of this process is given by 
Aji ■ iij, where yij denotes the number density of atoms in the ex- 
cited state j, and Ay, is the Einstein coeflicient for a spontaneous 
transition from level j to /. The total emissivity of the transition 
is 



e{T,ne) - hvji ■ Aj 



G(T) ■ 4 



(1) 



"■] 

where tig is the electron number density and G{T) the tempera- 
ture dependent contribution function 



GiT) = hvji ■ Aji 



llg ■ Tlj 



n 



nd_ nH_ 

rigi 11}] «(, 



(2) 



that comprises the relative abundance of the ionic species, the 
abundance of the element relative to hydrogen, and the number 
density of hydrogen relative to that of the electrons. These terms 
are accessib le thro ugh atomic data b ases, such as CHIANTI 
(Dere et al., 119971: Landi et al., I2006h . and can be easily eval- 
uated once the temperatures and densities are known. A rather 
strong assumption that has been used in th ese calculat i ons is 
to assume ionization equilibrium. However, iPeter et aP (l2006l) 
showed that ionization equilibrium is a fairly good approxima- 
tion that is valid for most parts in the simulation box. 

Nonetheless, it must be noted that future numerical experi- 
ments at higher resolution capable of resolving steeper temper- 
ature and velocity gradients will have to overcome this simplifi- 
cation. This has already been done in ID models, but is not yet 
feasible in 3D computations. 

The MHD calculations were performed using a nonequidis- 
tant grid in the vertical direction with a spacing of neighbor- 
ing grid points down to 75 km in the transition region, but the 
temperature resolution in these data is still insufficient for cal- 
culating the line emissivities in the transition region. To reach 
a temperature resolution below the typical contribution func- 
tion widths of transition region lines (log(r/[K]) w 0.2-0.3), the 
MHD results are interpolated in the vertical direction. The ver- 
tical distance between grid points in the transition region is thus 
reduced down to 40 km in the transition region and temperature 
gradients of 0. 1 in log T per grid cell are obtained on average 
at each height. All parameters discussed in the following were 
derived from such spatially interpolated data. 

The Doppler shift of the respective line (cf. Table [TJ is de- 
rived by calculating the first moment of the line profile. The av- 
erage line shift is obtained by averaging the line profiles hori- 
zontally before integrating along the line-of-sight. 

Table [T] shows the wavelengths and formation temperatures 
of the lines synthesized in this study. The lines are chosen so 
as to cover the entire transition region in temperature, ranging 
from T ^ 30.000 K to T ^ 1.3 x lO*" K. Also, these lines can 
be obser ved wi th instruments like SUMER/SOHO (e.g., Peter 
& Judge. [19991) . Temperatures labeled log Tii„e are the temper- 
atures, where the term ^T^'^^exp^-^), i.e., the product of 



Table 1. Emission lines synthesized in this study together with 
their wavelengths and line formation temperatures. 



line 


wavelength! 


line formation temperature 




[A] 


iog(r,„„/[K]) 


iog(r„„,„„/[K]) 


Sin*" 


1533.4 


4.44 


4.5 


Cn(" 


1334.5 


4.59 


4.6 


Si IV 


1393.8 


4.88 


4.9 


Cm 


977.0 


4.91 


4.8 


Civ 


1548.2 


5.02 


5.0 


Oiv 


1401.2 


5.19 


5.2 


Ov 


629.7 


5.40 


5.4 


Ovi 


1031.9 


5.48 


5.5 


Ne vm 


770.4 


5.80 


5.9 


Mgx 


624.9 


6.05 


6.1 



Notes. *" The formation temperatures for Si ii and C ii as listed here are 
most likely overestimated (see end of Sect. l3.lt 



maximum ionization fraction and excitation rate, of the respec- 
tive line has its maximum. The temperatures labeled log Tcomrib 
indicate the temperatures, where the maximum contribution is 
found in the simulation box (i.e., bulk of the emission in the re- 
spective line originates from regions with this temperature). The 
good correspondence between log Tu„e and log Tcomrib reflects 
that the transition region is basically isobaric in the simulation. 

Owing to the strong heating at the footpoints of the corona 
and the resulting strong increase of density towards the chromo- 
sphere, cool lines like Si ii and C ii tend to be formed well below 
their ioniza tion equilibr i um tem per ature. This issue is discussed 
in detail bv iJudee et aP (l2003l) and lPeter etaP (l2006l) . For con- 
sistency with the previous investigations outlined in this paper, 
e.g., when plotting Doppler shift vs. line formation temperature, 
we refer to the line formation temperatures listed in Table [Tlfor 
all ions, including Si ii and C ii. 

3.2. Doppler shifts derived from the MHD model 

Average Doppler shifts for each emission line of Table [T]are cal- 
culated from the emission line profiles at each time step. The line 
profile functions at each grid point are assumed to be Gaussian 
shaped and of the form 



where 



^/na>,h 



CO, I, 



exp 



(V - Vlos)' 



2kBT 



(3) 



(4) 



is the thermal width, I,,,, is the emissivity (Eq.[TJ at the respective 
grid point, and the vertical direction is assumed to be the line- 
of-sight. To obtain the Doppler shift of each emission line, the 
line profiles of all grid points are averaged horizontally at each 
height, and the average line profile is integrated along the line- 
of-sight. The first moment of this function corresponds to the 
shift of the average spectrum, i.e., the net Doppler shift of the 
respective emission line (cf. Fig.|2]i. 

The investigation of these average Doppler shifts shows 
that, in contrast to earlier work (e.g., Peter et al., l2004i |2006|) . 
a much better representation of the variation in average red- 
shifts with line formation temperature is obtained for the sim- 
ulation run of our model setup up to log(r/[K])«5.6, i.e., the 
low corona. Figure |2] shows a comparison between observed 
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Fig. 2. Average Doppler shifts derived from MHD model as 
function of line formation temperature. The height of the bars 
indicates the width of the distribution obtained by averaging hor- 
izontally over all line profiles at each height in the simulation 
box. 



Fig. 3. Temporal evolution of the relative change of mass in a 
given temperature interval as indicated in the legend. The av- 
erage mass over time in the respective temperature interval is 
denoted by < mass >. 



(data from Peter & Judge, Il999h and modeled Doppler shifts 
that show a close correlation at transition region temperatures 
up to log(r/[K])=5.6. Due to the high spatial resolution of the 
simulation, the dynamics in the simulation box are better re- 
solved, especially in the transition region where the temperature 
gradients are large. Also, the setup of the numerical experiment 
is more representative of the Sun outside of active regions (for 
which the transition region redsh ifts are well documented) than 
it was in earlier studies by, e.g., iGudiksen & Nordlundl (l2005h 
(see section 2). However, the observed blueshifts in the coronal 
part at log(r/[K])> 5.7 are not reproduced correctly by the MHD 
model owing to the impen etrable upper boun dary of the simula- 
tion box. As discussed by iPeter et aP (l2004h . this upper bound- 
ary condition is expected to have no influence on the transition 
region redshifts produced in the magnetically closed structures 
in the computational domain, but might prevent the formation of 
blueshifts in the coronal lines in the framework of our numerical 
experiment. On the real Sun, magnetically open regions might 
be connected to the solar wind and thus be responsible for the 
blueshifts of coronal lines. However, the upper boundary condi- 
tion of the model prevents the outflows and thus could suppress 
the observed blueshifts of coronal lines. 

4. Understanding transition region Doppler shifts 

The compliance between the observed transition region redshifts 
and the average Doppler shifts in the MHD model (cf. Fig. O 
allows the investigation of the origin of these line shifts based 
on our MHD model. 

4. 1 . Coronal mass balance 

As a first approach to tackling the line shifts in the model, 
the temporal variation in the total mass, i.e., the overall num- 
ber of particles, in various temperature intervals in the simula- 
tion box is determined. The temperature interval is chosen to be 
dlog{T/[K])-Q.l. Figure |3] shows examples for the variation in 
the mass in temperature intervals ranging from log(r/[K])=4.8- 
4.9 to log(r/[K])=5.6-5.7. The mass enclosed fluctuates between 
2% and 5% during the 30-minute simulation time. 

In principle, this small change in mass (or number of par- 
ticles A^) within a temperature interval could be due to an av- 



erage down- (or up-)flow, which could cause the observed av- 
erage downflow. However, an estimate of the average vertical 
velocity based on the observed variation in the number of par- 
ticles S = (N - {N})/{N} ~ 5%, the transition region scale 
height H^SMm, and the timescale of 15 min implies (v) = 
6-N/{nAT) - 6-H/t x 0.3 km/s, a value much too low to explain 
the observed redshiftsQ 

Because the intensity of optically thin emission lines scales 
with the density squared, the Doppler shifts can be approximated 
by {p^V;)/{p^), where (...) represents the average in a small 
temperature interval. This intensity-weighted velocity, of course, 
reflects the trend of the average Doppler shifts with line forma- 
tion temperature, as illustrated in Fig.|4] 

In Fig. m the velocities and densities that are evaluated are 
confined to regions in the simulation box within the respective 
temperature interval and above a height of z = 4.8 Mm. The 
minimum temperature at this height is log(r/[K])=3.8. It is as- 
sumed that plasma with temperatures below this threshold cor- 
responds to chromospheric material (cf. Fig. [1} rather than to 
transition region and coronal material. Since the purpose of this 
study is to investigate only the velocities of the transition re- 
gion and coronal plasma, chromospheric temperature areas are 
not considered. Otherwise, the material on average would be at 
rest, and a zero net mass flux would be observed. For tempera- 
tures above log(r/[K])=6.0, the small number of data points in 
the simulation box would result in unreliable predictions of the 
model, which is why such temperatures are not considered in the 
analysis. 

Cool material can reach high altitud es through pressur e im- 
balances or local dynamics, as shown bv lPeter et al.l (l2006h . e.g., 
in their Fig ly, S. These upwelling cool plasma fingers are not a 
good representation of spicules, because they lack the fine struc- 
ture and dynamics of observed spicules. These numerical exper- 
iments are not realistic enough to produce spicules, nonetheless 
they show the elevation of cool material into the corona. 

4.2. Pockets of cool gas 

A close inspection of Fig.|4]shows an upward directed intensity- 
weighted velocity at temperatures below log(r/[K])=4.0. More 

' N = nAH is the total number of particles in the corona, n = 10" 
m"^ is the particle density in the transition region, and A is a measure 
for the area corresponding to a horizonal cut through the box. 



4 



p. Zacharias et al.: Mass flows in the transition region 




logio line formation temperature [K] 



Fig. 4. Weighted velocity and mass flux as function of temper- 
ature in the simulation box. The solid line shows the intensity- 
weighted velocity (Sect. 14.1] ). This matches the Doppler shifts 
from the synthesized spectra (bars, for comparison, taken from 
Fig.|2]i. The average mass flux in the vertical direction is shown 
as a dashed line. 



important, the particle flux {n Vj) shows a change in sign at 
r^ilO'^K (dotted line in Fig.[4|. The upwards directed mass 
flux at low temperatures {T < 10"*K) is comparable in magni- 
tude to the downwards directed mass flux at temperatures above 
T - 10"* K. This would suggest that cool material is supplying 
the mass into the corona, where it is heated and then rains down. 
This is further supported by the finding that upflows (V; > 0) are 
identified up to heights of 10 Mm and 15 Mm for temperatures 
of log(r/[K]) < 4.0 and log(r/[K]) < 5.0 in the simulation box. 
These findings indicate a process by which cool pockets of mate- 
rial are pushed upwards from the chromosphere into the corona, 
while being heated slowly. The cool pockets are not visible at 
transition-region and coronal temperatures, since the heating oc- 
curs on a somewhat slower timescale than the rise of the mate- 
rial. Therefore, they cannot be observed as blueshifts in transi- 
tion region and coronal lines. This appears to be a phenomenon 
similar to the idea that cool spicular material is actually sup- 
plying mass to the corona t hat is raining down after it has been 
heated, as proposed, e.g., bv lAthav & Holzerl(ll982b . 

4.3. High-pressure plug of material that relaxes towards 
equilibrium by expelling material 

iHansteen et al.l ( 1201 Ol) suggest that the blueshifts in the low 
corona (log(r/[K]) > 5.6) can be explained as a natural con- 
sequence of episodically heated models, where the heating per 
particle is concentrated in the upper chromosphere, transition re- 
gion, or the lower corona. When a heating event occurs in the up- 
per chromosphere or transition region, cooler material is heated 
rapidly to coronal temperatures. In the lower corona, thermal 
conduction increases rapidly with rising coronal temperature, 
and in the case of a heating event, the chromospheric material 
below can be heated by a thermal wave. If these processes hap- 
pen expeditiously, the plasma will be heated more or less "in 
place" before it has time to move. Thus, newly heated transition 
material will be characterized by pressures higher than hydro- 
static, and the high-pressure material will expand into regions of 
lower pressure on either side: downwards towards the chromo- 
sphere and upwards towards the corona. 

We have found evidence for such a process in one of our 
simulation runs although on a larger scale. Plasma is ejected as a 



consequence of increased heating and high pressure in the lower 
transition region. Results are going to be presented in a second 
publication. This scenario as prop osed byHlansteen et al. (201^) 
or earlier bv lSpadaro et al.l (l2006l) is quite different from the one 
discussed above in Sect. 4.3. Probably, both processes occur in 
parallel, i.e., upwards accelerated cool parcels of plasma, as well 
as rapidly heated material in the low atmosphere. Future studies 
will have to show which, if any, of these two processes is respon- 
sible for the persistent Doppler shifts in the transition region. 

4.4. Field line related mechanism causing upflows and 
downflows 

So far, only vertical flows of material in the simulation domain 
have been investigated (Sect. 14. il l. The following discussion will 
focus on plasma flows along magnetic field lines. 

In principle, magnetic field lines can be traced easily in the 
simulation domain by following the magnetic field vector at a 
given grid point for a given distance. However, it is not straight- 
forward to follow such a magnetic field line over time, and some 
assumptions have to be made in order to develop a suitable 
technique. It is assumed that the field lines and the plasma are 
"frozen-in" due to the high electric conductivity of the plasma 
in the solar atmosphere. Thus, the velocity vector of the plasma 
at the corresponding time step and site (e.g., the top of the field 
lines) can be used to determine the displacement of the starting 
points of the field lines for the next time step. The highest point 
of each field line is used, because it is assumed that the field lines 
are more widely separated from each other the higher they are up 
in the atmosphere. Thus, errors resulting from jumping between 
different field lines can be minimized. The magnetic field vector 
is tracked in three dimensions from the starting point until the 
field line returns to the bottom or leaves the simulation box at 
the top. Field lines that leave the simulation box on either side 
can still be tracked as periodic boundary conditions are applied 
in the horizontal direction. The field line tracking is started in 
the regions of strong blueshift {\'o < -8 km/s) and in the layer 
of highest C iv (1548 A) emission in the simulation box. 

The finite electric conductivity in our numerical experiment 
means that working with "frozen-in" field lines is only an ap- 
proximation. The velocity of the slippage between magnetic 
field and plasma can be expressed by the term (///r)'^^, and is of 
the order of 1 km/s, calculated on the basis of the magnetic diffu- 
sivity T/ailO'" m^s"' and the time interval raiSOmin over which 
the field line is being followed. The calculated velocity of 1 km/s 
is below the velocity of the "frozen-in" field lines (perpendicu- 
lar to the magnetic field), and therefore the above assumption is 
justified. 

Field lines near a small X-type reconnection site have been 
selected for further analysis. The temporal evolution of a single 
field line is shown in a 3D view in Fig.|5]for four minutes (t=8- 
12 min) of simulation time. The vertical velocity v. along the 
magnetic field line is determined at each time step, since this is 
what the Doppler shift would be showing if the box is viewed 
from the top (situation at disk center). In Fig. |6] is shown at 
the location of highest C iv ( 1 548 A) intensity along the field line 
(since this is the location representing the observable Doppler 
shift). A quantitative analysis of the velocities shows a change 
in sign as the field line reconnects (from blue to red. Fig. |5)- 
In a first phase (t=0-8 min), a positive velocity is observed at 
the location of highest intensity along the line corresponding to 
an upflow of material. Whenever the field line reconnects, the 
vertical velocity switches signs. Plasma is draining out of the 
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Fig. 5. Temporal evolution of a single, selected magnetic field 
line near X-type reconnection site for four minutes simula- 
tion time. In addition, the vertical magnetic field component B; 
scaled from -1000 G (black) to -i-lOOO G (white) is shown at 
height zero (z-0). 



loop during the period of time when the field line is connected 
to a region different than in the first phase. Another connectivity 
change in the field line occuring in the time period between t=12 
min and t=13 min is not shown in Fig.|5] The vertical velocity 
switches signs to positive again (cf. Fig. |6l), corresponding to 
an upflow of material within the loop. Vertical velocities of two 
spatially fixed regions of mass loading (v^ > at t=6 min) and 
draining (v; < at t=10 min) are indicated by the red and blue 
dotted lines in Fig.|6] The temporal evolution of these velocities 
shows Uttle variation, indicating that these sites remain rather 
stable during the simulation. 

Figure [T] shows a schematic view of the field line configu- 
ration during the upflow (indicated by 1, blue) and downflow 
(indicated by 2, red) phases. The phases labeled 1 and 2 cor- 
respond to phases named identically in Fig. |6] and are used to 
illustrate the two-step process outlined above. The coronal part 
of the magnetic field line is initially connected to a region below, 
where the material is being heated and accelerated. This leads to 
an upflow of material into the corona that is being observed as a 
blueshift along the loop (phase 1). After magnetic reconnection 
has led to a reconfiguration of the field, the magnetic field line 
is connected to a different chromospheric region with less heat- 
ing from below. As a consequence, the material flows back down 
along the loop leading to the observed redshift (phase 2). A net 
redshift is observed, since the upflows and downflows take place 
at different temperatures and densities (cf. section l4Tt . 

5. Discussion and conclusions 

This study of the origin of Doppler shifts in a 3D MHD model 
extending from the photosphere to the corona reveals persistent 
redshifts at transition region temperatures. The comparison of 
the modeled to the observed Doppler shifts shows a close simi- 
larity between both quantities with large apparent downflows at 
transition region temperatures. The relation between the Doppler 
shifts and the mass flux in the model can be analyzed as a next 
step. The synthesized Doppler shifts are reproduced well by an 
"intensity weighted" velocity, i.e., the velocity weighted by the 




time [min] 



Fig. 6. Temporal evolution of the vertical velocity v. along the 
field line followed in Fig.|5]at the position of highest C iv (1548 
A) intensity. (1) and (2) indicate the two different field line con- 
figurations shown in Fig. |5] associated with different flow pat- 
terns along the field line. (1): Material is flowing upwards along 
the field line at the site of highest C iv (1548 A) corresponding to 
an observed blueshift. (2): Material is flowing downwards along 
the field line at the site of highest C iv intensity corresponding to 
an observed redshift. 




Fig. 7. Two-dimensional illustration of the changing magnetic 
field line configuration shown in Fig. |5] The phases labeled 1 
and 2 correspond to phases named identicaUy in Fig.|6] 



density squared, {p- v)/{p-). The basis for this is that the emis- 
sivity is proportional to the density squared (oc p^). 

The coronal mass is approximately constant in these simu- 
lations, varying between 0.2% and 0.5% on average. This varia- 
tion only cannot explain the persistent transition region Doppler 
shifts that are a few km/s. This discrepancy is mitigated by the 
finding that a simple interpretation of the Doppler shifts in terms 
of an average vertical mass flux leads to a systematic overesti- 
mation of the mass flux, an effect that is most obvious in the 
upper transition region. The presence of upflows at temperatures 
below log(r/[K])=4.0 indicates a process by which cool pock- 
ets of material are pushed upwards into the corona from below. 
As they rise, they heat up slowly. Due to their low temperature, 
they are not observable as blueshifts in the transition region and 
coronal lines. 

When investigating the magnetic field topology, there is ev- 
idence for a scenario in which magnetic field lines are inter- 
mittently connected to subjacent regions where heating takes 
place on different scales. This leads to either upflows observed 
as blueshifts or downflows observed as redshifts along the field 
lines. The regions of mass loading and draining are fairly stable 
in space and time and can thus serve as mass reservoirs or sinks. 
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Plasma on field lines traversing these regions experiences either 
an upflow or a downflow depending on the connectivity of the 
field line. 

As a conclusion, the good correlation between observed 
and model-synthesized transition region Doppler shifts below 
log(r/[K])=5.6 is due to cool pockets of plasma propelled up- 
wards and heated to transition region and coronal temperatures 
after wards, a process s imilar to the one previously suggested by, 
e.g., lAthav & Holzeil (fT98l . This is in contrast to recent work 
by Hansteen et al. (2010), who suggest that chromospheric ma- 
terial is heated fast in the low atmosphere and expands into the 
corona afterwards. Future investigations will have to decide on 
the merits of these alternatives. 

Acknowledgements. We thank the referee, Phil Judge, for very useful com- 
ments on the manuscript. This research has been supported by the Deutsche 
Forschungsgemeinschaft (DFG). 



References 

Athay, R. G. 1984, ApJ, 287, 412 

Athay, R. G. & Holzer, T. E. 1982, ApJ, 255, 743 

Bingert, S. & Peter, H. 201 1. A&A, in press, arXiv 1 103.6042 

Brandenburg, A. & Dobler, W. 2002, CoPhC, 147, 471 

Cook, J. W., Cheng, C.-C, Jacobs, V. L., & Antiochos, S. K. 1989, ApJ, 338, 
1176 

Dere, K. P., Landi, E., Mason, H. E., Monsignori Fossi, B. C, & Young, P. R. 

1997, A&AS, 125, 149 
Doschek, G. A., Bohlin, J. D., & Feldman, U. 1976, ApJ, 205, L177 
Gudiksen, B. & Nordlund, A. 2005, ApJ, 618, 1020 
Hansteen, V. H. 1993, ApJ, 402, 741 

Hansteen, V. H., Hara, H., Pontieu, B. D., & Carlsson, M. 2010, ApJ, 718, 1070 
Hassler, D. M., Rottman, G. J., & Orrall, E Q. 1991, AdSpR, 11, 141 
Judge, R G., Carlsson, M., & Stein, R. F 2003, The Astrophysical Journal, 597, 
1158 

Landi, E., Del Zanna, G., Young, P R., et al. 2006, ApJS, 162, 261 

Maiiska, J. T. 1992, The Solar Transition Region (Canbridge University Press) 

Peter, H. 1999, ApJ, 516, 490 

Peter, H., Gudiksen, B., & Nordlund, A. 2006, ApJ, 638, 1086 
Peter, H., Gudiksen, B. V., & Nordlund, A. 2004, ApJ, 617, L85 
Peter, H. & Judge, P G. 1999, ApJ, 522, 1 148 

Spadaro, D., A. F Lanza, J. T. K., & Antiochos, S. K. 2006, ApJ, 642, 579 
Spitzer, L. 1962, Physics of fully ionized gases (Interscience, New York) 
Wilhelm, K., Curdt, W., Marsch, E., & et al. 1995, Sol. Phys., 162, 189 



